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The observational indications for a constant specific star-formation rate (sSFR) 
in the redshift range z = 2 — 7 are puzzhng in the context of current galaxy- formation 
models. Despite the tentative nature of the data, their marked conflict with theory 
motivates a study of the possible implications. The plateau at sSFR~2 Gyr""'^ is hard 
to reproduce because (a) its level is low compared to the cosmological specific accretion 
rate at 2 > 6, (b) it is higher than the latter at z ^ 2, (c) the natural correlation between 
SFR and stellar mass makes it difficult to manipulate their ratio, and (d) a low SFR 
at high z makes it hard to produce enough massive galaxies by z^2. Using a flexible 
semi-analytic model, we explore ad-hoc modifications to the standard physical recipes 
trying to obey the puzzling observational constraints. Successful models involve non- 
trivial modifications, such as (a) a suppressed SFR at 2; > 4 in galaxies of all masses, 
by enhanced feedback or reduced SFR efficiency, following an initial active phase at 
z > 7, (b) a delayed gas consumption into stars, allowing the gas that was prohibited 
from forming stars or ejected at high z to form stars later in more massive galaxies, 
and (c) enhanced growth of massive galaxies, in terms of either faster assembly or 
more efficient starbursts in mergers, or by efficient star formation in massive haloes. 
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1 INTRODUCTION 

The bulk of the stellar mass observed today in galaxies is 
built up at high redshift, where star formation and mass as- 
sembly are very efficient. This makes observations at high 
redshift extremely important in understanding galaxy evo- 
lution in general. For example, the maximum value of the 
mean cosmological star-formation rate (SFR) density is 
achieved at redshift 1-3 (e.g. Hopkins & Beacom 2006). 

Although high redshift observations are naturally 
plagued with larger uncertainties than local data, it is clear 
that a successful model of galaxy formation and evolution 
should match them within realistic error margins. Current 
models are often mainly tuned to reproduce the low redshift 
universe, and it is therefore very instructive to test them 
against high redshift observations. For example, Fontanot 
et al. (2009) and Guo et al. (2010) point out that the low 
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mass end of the stellar mass function is built too quickly at 
high redshift in current semi-analytic models (SAMs here- 
after), while Khochfar et al. (2007) claim that their SAM 
can reproduce the evolution of the faint end of the lumi- 
nosity function. It seems however likely that the evolution 
of the stellar mass function over time alone is insufficient 
to fully constrain the models, as was argued by Neistein & 
Weinmann (2010). Consequently, it is a crucial next step 
to compare the observed SFR of high redshift galaxies to 
models in more detail than previously done, where the main 
focus was in trying to match the global star formation rate 
density. 

Pioneering observational estimates of the SFR (rhstar) 
and stellar mass (mstar) indicate that the specific SFR 
(sSFR, /mstar) is roughly constant in time throughout 
the redshift range 2: = 2 — 7, for galaxies with roughly the 
same mass matar ~ (O.2-1)-1O^°M0, at alevel sSFR ~ 1-2 
Gyr~^. (e.g. Stark et al. 2009a; Gonzalez et al. 2010; Labbe 
et al. 2010a, 2010b). There are indications that this sSFR 
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plateau is associated with a rather constant sSFR within 
each galaxy as it grows (Papovich et al. 2010; Stark et al. 
2009a). The plateau is hard to reconcile with the current 
theoretical wisdom for several reasons. 

In current models of galaxy formation the high-z SFR 
is assumed to a large extent to be driven by the fresh gas 
supply (see Keres et al. 2005; Dekel et al. 2009; Bouche et 
al. 2010; and basically all the SAM; but see also Narayanan 
et al. 2010 who discuss a potential merger origin for the 
population of submillimeter galaxies at 2:~ 2). The specific 
cosmological accretion rate of baryons is steeply declining 
with time, M/M oc {l + zf-^= (Neistein & Dekel 2008; Dekel 
et al. 2009). In particular, at 2 ~ 7 the specific accretion rate 
is higher than the observed sSFR by a factor of a few, and 
at « ~ 2 it is lower than the sSFR by a similar factor. This 
is in marked contrast to the observed plateau in the sSFR. 
A constant sSFR during the evolution of each galaxy (main 
progenitor) would require either an exponential growth in 
time of TTistar and mstar (this is if most stars are formed in 
situ to the main progenitor), or a non-trivial combination of 
effective SFR and stellar assembly rate as a function of time 
and mass. An obvious related difficulty is introduced by the 
fact that a low SFR at high z could make it hard to produce 
enough massive galaxies by 2 ~ 2 to match the bright end of 
the observed galaxy mass function at that epoch. 

Given the marked contrast between the indicated obser- 
vation and the current models of galaxy formation, we ap- 
peal to a special semi-analytic tool. Traditional SAMs (e.g 
Kauffmann et al. 1993; Cole et al. 2000; De Lucia & Blaizot 
2007) describe the processes that are responsible for galaxy 
evolution by physically motivated recipes that are fixed a 
priori. They are thus geared to solve the 'forward problem', 
i.e., test to what extent the assumed set of physical recipes 
provides a match to the observed properties of the galaxy 
population. This methodology is not ideal for exploring a 
large variety of physical recipes, some of which may need 
to deviate significantly from the standard assumptions. Our 
approach here is to solve the 'inverse problem', where we 
investigate how the observations at high z constrain the ba- 
sic recipes, especially those associated with the processes of 
mergers, star formation, and feedback. This is achievable 
using the method of Neistein & Weinmann (2010, NWIO 
hereafter), which allows freedom in choosing the recipes of 
interest, and a very efficient exploration of a broad param- 
eter space. 

In order to reproduce the sSFR plateau, we will try to 
either lower the SFR efficiency after a very short period of 
high efficiency at z > 7, or to enhance the suppression of 
SFR by feedback at z > 4, especially in high-mass haloes, 
to be followed by an enhanced SFR in the retained or rein- 
corporated gas at 2: ~ 2 — 3. In both cases, the low sSFR at 
high z makes it difficult to form enough massive galaxies at 
2 ~ 1 — 3, unless the rate of mass assembly due to mergers 
and the associated starbursts are pushed to their limits. 

The models presented in this paper are not a priori 
physically motivated. However, we have tried to keep them 
simple, and to minimize deviations both from the standard 
model and from a monotonic dependence on halo mass and 
time. We also try to keep the models as physically plausible 
as possible. Our aim in this paper is not to find the "right" 
model. The main goal is to investigate the approximate na- 



ture of the needed changes, and how they impact on the 
properties of the galaxy population other than the sSFR. 

The measurements of SFR and stellar mass at high 2 
are still at their infancy, and therefore their interpretation, 
in particular the sSFR plateau, is uncertain and highly con- 
troversial. The main source of uncertainty is the obscuration 
by dust, where several authors bring convincing arguments 
for little or no dust at 2 > 4 (e.g. Bouwens et al. 2009; 
Finkelstein et al. 2010), while others do apply dust correc- 
tions and obtain higher values of sSFR. Stellar masses may 
also be affected by systematic errors. For example, uncer- 
tainties in the treatment of the TP AGB-phase may lead to 
errors in the stellar mass estimates at 2 = 2 — 3 by factors 
of ~ 2 (e.g. Maraston et al. 2006; Magdis et al. 2010a). The 
disagreement between a constant sSFR and our current the- 
oretical wisdom is so pronounced that it motivates a study 
of the theoretical implications despite the observational un- 
certainties, adopting the validity of the sSFR plateau as a 
working assumption. An alternative way out from the puz- 
zle might be to assume a time-dependent stellar initial mass 
function (IMF) (e.g. Dave 2010), but this is kept beyond the 
scope of the current paper. 

The outline of the paper is as follows. In section [2] we 
summarize the observed sSFR at high redshift and explain 
the points of tension with theory. In section [31 we explain 
the NWIO method that we use here. In section |4l the main 
part of this paper, we demonstrate how a simple standard 
SAM fails to reproduce the sSFR plateau, and proceed by 
making controlled changes to this model in order to bet- 
ter reproduce the observations. In sections 15.11 and 15.21 we 
discuss possible physical motivations for these changes. In 
section 15.31 we comment on the observational uncertainties 
regarding obscuration by dust. In section lS^ we compare our 
results to previous work. Finally, in section[Sl we present our 
conclusions. Throughout the paper, we refer to the redshift 
range 2 = 2 — 3 as "intermediate" redshift, to 2 = 4 — 6 as 
"high" redshift, and to 2 > 6 as "very high" redshift. Our 
models are based on dark-matter merger trees from the Mil- 
lennium Ai'-body simulation (Springel et al. 2005), and we 
thus use WMAPl cosmological parameters. In particular, 
we quote masses assuming h — 0.73. 



2 THE PROBLEM 
2.1 Observations 

Acorrelation between stellar mass and SFR has been ob- 
served at various redshifts up to 2 ~ 4 (e.g. Daddi et al. 
2007; Stark et al. 2009; Karim et al. 2011). Figure [T] shows 
a compilation of observational estimates of sSFR as a func- 
tion of redshift for star-forming galaxiej^ of a similar stellar 
mass ~ (0.2— 1) • 10^" Mq. The data reveal a rather constant 
sSFR ~ 2 Gyr^^ in the redshift range 2 = 2 — 7 (Feulner et 
al. 2005; Yan et al. 2006; Eyles et al. 2007; Daddi et al. 2007; 
Stark et al. 2009; Gonzalez et al. 2010; Labbe et al. 2010a,b; 

^ typically selected at high redshift as Lyman-break galaxies 
(LBGs, Steidel et al. 1999), which possibly excludes a popula- 
tion of low SFR galaxies (e.g. Richards et al. 2011). This will 
however only increase differences between model and basic theo- 
retical predictions. 
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Figure 1. The observed sSFR plateau. Shown are measurements 
of specific star-formation rate as a function of redshift for galaxies 
in a similar stellar mass range ~ (O.2-l)-lOl°M0. The references 
are marked and listed in the text. The grey belt captures most 
of the measurements, and reflects the uncertainty of ±0.3 dex 
estimated by Gonzalez et al. (2010). It indicates a constant sSFR 
~ 2 Gyr~^ in the range z = 2 — 7, followed by a steep decline 
toward 2 = 0. All references given in square brackets refer to a 
rather high median mass of ~ 10^° Mq. 



Rodighiero et al. 2010; Gonzalez et al. 2010; IVIagdis et al. 
2010a,b; McLure et al. 2011). except for three higher esti- 
mates (Yabe et al. 2009; Schaerer & de Barros 2010, Shim 
et al. 201lj3- The main reason for these higher estimates is 
the larger correction for dust extinction assumed by these 
authors (to be discussed in section [53)) . as well as different 
treatments of nebular emission lines and different assump- 
tions concerning star-formation histories (to be discussed in 
section r4.6|) . Additionally, Shim et al. (2011) only include 
galaxies with indications for Ha emission, which will bias 
the estimate of the sSFR high. We note that all the esti- 
mates above do not include submillimeter galajcies, which 
are outliers to the relation between stellar mass and SFR, 
simply because these tend to have stellar masses above the 
limit we consider here (e.g. Daddi et al. 2007). At 2 < 2, 
the sSFR declines steeply (Noeske et al. 2007; Elbaz et al. 
2007; Dunne et al. 2009; Oliver et al. 2010; Rodighiero et 
al. 2010). The grey belt in Fig.[T]tries to capture the overall 
trend, reflecting an uncertainty of ±0.3 dex as estimated by 
Gonzalez et al. (2010), and ignoring the two high estimates. 
As will be discussed below, this observed sSFR plateau is 
puzzling — its level is surprisingly high at z ~ 2 and sur- 
prisingly low at z > 4. For the purpose of the theoretical 
analysis of the current paper, we adopt the sSFR plateau as 
marked by the grey belt. 



^ We obtained part of the estimates by dividing the median SFR 
by the median stellar mass (for Yan et al. 2006; Eyles et al. 2007; 
Yabe et al. 2009). In some other cases, the estimates are based 
on an extrapolation of the sSFR-stellar mass relation to a stellar 
mass of ~ O.5-lO"'Af0 (Daddi et al. 2007; Rodighiero et al. 2010). 



Figure 2. Evolution of sSFR in the SAMs of NWIO, and of 
the specific dark matter accretion rate. Shown are four different 
models (I, II, III and V, with results of models II and V being 
indistinguishable and thus represented by one line) by the SAM 
of NWIO for galaxies in the mass range [2-10^, 10^''] Mq (curves 
in colour). The completeness limits trying to mimic the observed 
ones are described in the text. In all models the sSFR is steeply 
declining in time, not reproducing the observed sSFR plateau 
marked by the grey belt from Fig. 1. Also shown is the specific 
dark matter accretion rate onto haloes of log(Mijaio) ~ IO^^Mq 
according to Neistein & Dekel (2008) (dashed black line), and the 
same quantity multiplied by a factor of 2, to account for the effect 
of instantaneous mass loss from newly formed stars, as assumed 
in the models (solid black line). 

2.2 Tension with theory 

Here, we outline the main potential points of tension be- 
tween the observed sSFR and theoretical predictions both 
from relatively detailed SAMs and simple analytical argu- 
ments. 

2.2.1 Tension with SAMs 

In Fig. O we show the sSFR as a function of redshift at 
a fixed mass for four of the models presented in NWIO, in 
comparison with the observed sSFR plateau. To account for 
the observational completeness limit as indicated by Stark 
et al. (2009), we only take into account model galaxies with 
log(SFR) > 0.25, 0.45 and 0.5 yr"^ at z=4, 5 and 6 respec- 
tively. The four models are described in detail in NWIO with 
the same numbering used here and can be summarized as fol- 
lows: I) model without SN feedback, II) model without ejec- 
tive SN feedback, HI) model including cold accretion [which 
is the model most similar to other standard SAMs], and V) 
model in which cooling and star formation shuts down af- 
ter major mergers. All of those models have been tuned to 
reproduce key observables like the stellar mass functions at 
different redshifts, and star formation rate at z=0. Remark- 
ably, all models show an extremely similar behaviour despite 
their fundamental differences, all in disagreement with ob- 
servations. They overestimate the sSFR by about an order 
of magnitude at z ~ 6, and underestimate it by around 0.3 
dex at 2 ~ 2. 

Other current SAMs show a similar behaviour. For ex- 
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Figure 3. Evolution of sSFR in variations of the standard SAM 
used in this work, and of the specific dark matter accretion 
rate. The predicted dark matter accretion rate onto haloes with 
log(Mhaio) ~ 10^'^ Mq according to Neistein & Dekel (2008) 
(dashed black line) and multiplied by a factor of 2 (black line), 
compared to the standard model (blue dashed line, described in 
section r4.1l l. the standard model without feedback (red solid line), 
and the standard model without feedback and stellar mass loss 
(green dot-dashed line). The sSFR in the raw model without feed- 
back and mass loss matches the total specific accretion rate. The 
mass loss adds a factor of two. 



for the galaxies considered here, this agreement is notewor- 
thy. We see in the figure that the feedback as implemented 
in the standard model does not have a significant effect on 
the sSFR at 2: > 4, while it gradually reduces the sSFR at 
lower redshifts. Compared to the models of NWIO shown in 
Fig. [21 the standard model that we use in this paper has a 
low sSFR ai z < 2, due to efflcient feedback at late times. 
This is not relevant for studying the plateau. 

fn what follows, we will demonstrate that despite this 
serious tension with theory, the sSFR plateau can in prin- 
ciple be reproduced by models of galaxy evolution, but it 
takes non-negligible modifications to common ingredients of 
these models. 



3 THE FORMALISM 

In this section we describe the formalism we use for mod- 
eling the evolution of galaxies. For more details about the 
methodology the reader is referred to NWIO. It is shown 
there that the results of our model are very similar to those 
given by a standard SAM, although the recipes are simpli- 
fied and schematic. In the context of this work, the simphc- 
ity of the model allows us to tune it easily, without losing 
the complex interplay between different process like dark- 
matter growth, cooling, SF, feedback, and merging. The 
code is available for public usage through the Internet (see 
http : //www.mpa-garching.mpg . de/galf orm/ sesam) 



ample, Lacey et al. (2010) show that the sSFR of galaxies 
in the Baugh et al. (2005) model have sSFR > 10 Gyr"^ 
at z—6, an order of magnitude higher than observational re- 
sults. Daddi et al. (2007) indicate that their observed SFR 
at z ~ 2 is significantly higher than the values predicted by 
the model of Kitzbichler & White (2007). Finally, we have 
confirmed ourselves that the model of De Lucia & Blaizot 
(2007) predicts results very similar to our model predictions 
shown in Fig. 2 (see also Guo & White 2008, their Fig. 3). 

2.2.2 Tension with basic theoretical considerations 

The observed sSFR plateau is in disagreement with the stan- 
dard wisdom concerning galaxy evolution. First, the average 
specific accretion rate into dark-matter haloes of a given 
mass is rapidly increasing with redshift, roughly in propor- 
tion to (1 -I- zf-^ (Neistein & Dekel 2008). Second, galaxies 
are more dense and gas rich at high-redshift, which is ex- 
pected to lead to higher SFR (e.g. Button et al. 2010). In 
agreement with these studies, Bouche et al. (2010) find in 
their idealized model that the sSFR of individual galaxies is 
indeed monotonically decreasing with time. To illustrate the 
first point, we show in Fig. [2]the approximation for the aver- 
age specific dark matter accretion rate M /M for haloes with 
log(Mhaio) ~ W^^Mq (Neistein & Dekel 2008), and compare 
it to the standard model, described in section \4A\ in Fig. [S] 
Remarkably, when we remove feedback from the standard 
model, it predicts a sSFR evolution in excellent agreement 
with twice the specific dark matter accretion rate. As shown 
in the figure, the remaining factor of two difference is fully 
explained by instantaneous stellar mass loss. Given that the 
ratio between stellar mass and halo mass is only about 5% 



3.1 Merger trees 

We use merger trees extracted from the Millennium A'^- 
body simulation (Springel et al. 2005). This simulation was 
run using the cosmological parameters (f2m, f^A, h, as) = 
(0.25, 0.75, 0.73, 0.9), with a particle mass of 8.6- 10* /i"^Mq 
and a box size of 500 h~^Mpc. The merger trees used here 
are based on subhaloes identified using the SUBFIND algo- 
rithm (Springel et al. 2001). They are defined as the bound 
density peaks inside fof groups (Davis et al. 1985). More 
details on the simulation and the subhalo merger-trees can 
be found in Springel et al. (2005) and Groton et al. (2006). 
The mass of each subhalo (referred to as Mh in what fol- 
lows) is determined according to the number of particles it 
contains. Within each FOF group the most massive subhalo 
is termed the central subhalo of this group. Throughout this 
paper we will use the term 'haloes' for both subhaloes and 
the central (sub)halo of FOF groups. 

3.2 Quiescent evolution 

Each galaxy is modeled by a 4-component vector, 

/ mstar \ 

m = , (1) 

mm 

\ mth / 

where mstar is the mass of stars, rricoid is the mass of cold 
gas within the disk, mm is the mass within cold filaments 
streaming within the host halo into the central galaxy, and 
jTifb is the mass currently made unavailable for star forma- 
tion by stellar feedback. We use the term 'quiescent evo- 
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lution' to mark all the evolutionary processes of a galaxy, 
except those related to mergers. 

All the models in this work assume that fresh gas is 
added to a galaxy only by cold filaments, increasing the 
mass of mm- The infall rate into filaments is assumed to be 
proportional to the dark-matter growth rate. 



Maceration = 0-17 Mh 



(2) 



Here 0.17 is the cosmic baryonic fraction, and Mh is the 

rate of dark-matter smooth accretion which does not include 
mergers with resolved progenitors (if Mh < we use a gas 
accretion rate of zero). 

The mass of cold gas within the disk is increased due to 
the free infall of cold filaments from the outer parts of the 
host halo. Wc mimic this effect by assuming that gas joins 
the disk with a specific rate /c. 



T-cold]ff — — ["lfil]ff — /c ■ fifil ■ 



(3) 



The efficiency /c = /c(Mh,t) is a function of the host halo 
mass Mh and the cosmic time t only, and is given in units 
of Gyr-i. 

We assume that the SF rate is proportional to the 
amount of cold gas. 



[rhstar]sF = ~ ['Tlcold]sF = fe ' Tlcold , 



(4) 



where /s = /s(Mh,t) is a function of the halo mass and 
time, in units of Gyr~^. For each SF episode we assume 
that a constant fraction of the mass is returned back to the 

cold gas component due to SN events and stellar winds. This 
recycling is assumed to be instantaneous, and contributes 



nicold 



recycling 



J recycling 



(5) 



Following NWIO, we use R — 0.5 for all models. This is 
the recycled fraction for a Chabrier (2003) IMF at 13.5 Gyr 
after a star burst according to the Bruzual & Chariot (2003) 
stellar population models. We note that this is the only point 
where the assumption on the IMF enters our model. 

Cold gas can be affected by feedback, which means that 
it becomes unavailable for star formation and moves from 
the cold phase to the feedback phase. Assuming that stellar 
feedback immediately follows star formation, this feedback 
should be in proportion to the SF rate, 



[mfbjfocdback 
fd [msta 



' [m,coldJfgejjback ' 
= /d/sWcold . 



(6) 



We model feedback by a function of halo mass and time, 
/d = /d(Mh,t). 

Once gas has been made unavailable for star formation 
due to feedback, we allow it to return to the cold phase, with 
a re-incorporation efficiency: 



[rHcoidli-c = - ['Tifb]„ = /rcmfb- 



(7) 



Note that the feedback mechanism we use is moving gas 
from the cold phase to the feedback phase. The mass 
within filaments is not participating in the feedback and re- 
incorporation processes. If the cold gas mass becomes neg- 
ative in a given timestep due to strong feedback, the star 
formation rates are adjusted such that a cold gas mass of 
zero is produced. 

To conclude, each process is described by one function 



which depends on the host halo mass and time only. All pro- 
cesses discussed in this section can be written in a compact 
form by using the following differential equations: 



m = Am -I- BMh , 



where 



/ 



V 



(1 " R)f. 
(1 - R)f. - /d/s 


/d/s 





/c 
-fc 







-/rc 



f » ^ 

0.17 

V / 



(8) 



(9) 



(10) 



Photoionization heating of the intergalactic medium is 
assumed to suppress the amount of cold gas available for SF 
within low mass haloes. This effect is critical for modeling 
the formation of dwarf galaxies. The minimum halo mass 
of ~ 2 ■ lO^'Vi^^M© in the Millennium simulation, which 
we use here, does however not allow a detailed modeling 
of small mass galajcies. Thus, instead of implementing a de- 
tailed treamtment of reionization, we simply assume that all 
the gas is kept hot until redshift 9, where cooling and SF 
are allowed to start. This is a higher redshift than in NWIO, 
which we found is needed to produce a high enough number 
of galaxies at very high redshifts. 

3.3 Mergers and satellite galeixies 

Satellite galaxies are defined as all galaxies inside a 
FOF group except the main galaxy inside the central (most 
massive) subhalo. Once the subhalo corresponding to a given 
galaxy cannot be resolved anymore, it is considered as hav- 
ing merged with the central halo. Due to the effect of dynam- 
ical friction, the galaxy is then assumed to spiral towards the 
center of the FOF group and merge with the galaxy in the 
central halo after a significant delay time. 

At the last time the dark matter subhalo of a satellite 
galaxy is resolved we compute its distance from the central 
halo (rsat), and estimate the dynamical friction time using 
the formula of Binney (1987), 



tdt = Qdf ■ 



1.17Kr,i 



Gmsat In (1 -I- Mh/rrisat) 



(11) 



For jTisat we use the baryonic (stars -|- cold gas) mass of the 
satellite galaxy plus the minimum subhalo mass which can 
be resolved by the Millennium simulation. Vy, Mh are the 
virial velocity and mass of the central subhalo. If a satellite 
falls into a larger halo together with its central galaxy we 
update tdi for both objects according to the new central 
galaxy. 

While satellite galaxies move within their FOF group, 
they suffer from loss of their extended gas reservoir due to 
tidal stripping. We assume that all satellite galaxies are los- 
ing their reservoir of filament gas exponentially, on a time 
scale of a few Gyr. In order to properly model this stripping 
we modify A by subtracting a constant an from one of its 
elements; 



A,at(3,3) =-/c 



ah ■ 



(12) 
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Figure 4. The cooling efRciency as a function of halo mass and 
redshift, as used in all the models described in sec. (4] 



Note that a constant in the diagonal of A gives an expo- 
nential time dependence. However, the actual dependence 
of JTifii on time for satellite galaxies is more complicated due 
to contributions from accretion. In general the parameter ah 
should depend on the dynamical time of the host halo. For 
simplicity we consider it to be a constant here. We assume 
that the gas which is in the feedback phase is not stripped. 

When galaxies finally merge we assume that a SF burst 
is triggered. We follow Mihos et al. (1994), Somerville et al. 
(2001) and Cox et al. (2008) and model the amount of stars 
produced by 



ArUa 



/burst (Wl, cold + "l2,cold) 



where 



/bur 



at 



mi 
m2 



(13) 



(14) 



Here rrii are the baryonic masses of the progenitor 
galaxies (cold gas plus stars), cold is their cold gas mass, 
and ai,, Oc are constants. 

The burst duration has been shown to vary in hydro- 
dynamical simulations between tens of Myr to a few Gyr 
depending on the merger mass ratio, and whether multi- 
ple bursts are considered or just the main peak (Cox et al. 
2008). We use a timescale of 10 Myr in all our models, fol- 
lowing De Lucia & Blaizot (2007). Merger-induced bursts 
cause feedback in the same way as quiescent star formation. 



4 MODIFYING A SIMPLE MODEL 

Below, we present a simple standard SAM, which we use 
as a starting point for our tuning procedure. This simple 
model includes features common to many current models of 
galax;y formation. It is kept as simple as possible in order to 
facilitate tuning and to simplify interpretation of changes to 
the model. In section 14.21 we present 6 alternative models 
which reproduce the sSFR plateau. 
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Figure 5. Star formation efRciency, Cs, as a function of redshift 
and halo mass, in the different models as indicated. Top: in the 
standard model and in models FBI and FB2, es is constant at 
all times and halo masses. In models SFl, SF2, and SF3, it is a 
function of redshift only. Bottom: in model SF4, is a function 
of both redshift and halo mass. 



4.1 The standard model 

The ingredients of this model were chosen based on the sim- 
plicity of the physical processes involved and consistency 
with observational constraints. 

• We assume that all accretion is cold, with /c — 1/tdyn, 
and that accretion is quenched at z < 3 above A/h ~ 1-6 • 
lO^^M© (e.g. Dekel & Birnboim 2006; Cattaneo et al. 2006; 
Ocvirk et al. 2008; Dekel et al. 2009). To avoid sharp breaks 
in the stellar mass function, we smoothed the transition both 
in redshift and in mass by hand, as shown in Fig. [J] 

• The rate at which cold gas is turned into stars is 

/s = Es/tdyn , (15) 

where es is the star-formation efficiency. For fdyn.disk, we 
assume 

3A • 7?vir 



^dyn.disk — 



\/2-Ki 



' 0.0072 • tHubbic , 



(16) 



where the halo spin parameter A = 0.03 (according to the 
mean value found by Munoz-Cuartas et al. 2010 in N-body 
simulations). We use a constant star formation efficiency 
Es ~ 0.01, as shown in Fig.[5]as blue line. This is comparable 
to the estimates of Krumholz & Tan (2007), and similar to 
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Figure 6. Feedback as a function of halo mass and redshift in our models. Top left: feedback efficiency in the standard model and 
all SF models. Top right: in model FBI. Bottom left: in model FB2. Contours mark log(/(j) = -0.5, 0, 0.5, 1, 1.5 and 2. Bottom 
right: reincorporation efficiency /rc in units of 1/tdyn in model FB2. The main new feature in the FB models is strong feedback at ^ > 4, 
and especially so in massive galaxies. Reincorporation is effective at z < 3 and in massive galaxies. 



the efficiencies used in other recent models (e.g. Krumholz 
& Dekel 2010; Bouche et al. 2010 and Guo et al. 2010). 

• We include a simple prescription for stellar feedback 
which is the usual way to reproduce the low mass end of the 
stellar mass function (Dekel & Silk 1986). The formation 
of stars aflects the surrounding cold gas, making part of 
it unavailable for star formation for a certain amount of 
time. The rate at which gas is made unavailable for star 
formation is equal to the SFR times an efficiency /d. As in 
most other SAMs, /d is a combination of a constant term 
(which is called the "reheating" part of feedback in De Lucia 
& Blaizot 2007), and a term that is inversely proportional 
to the virial velocity to some power, which should roughly 
mimic ejection of material due to SN explosions (e.g. Bower 
et al. 2006; De Lucia & Blaizot 2007; Guo et al. 2010). The 
feedback efficiency is thus given by 

/d=<5 + (^)-". (17) 

7 

We use a=3.5, 7=161 km/s, 5 = We do not let log(/d) 

^ A non-zero delta is necessary in order to have some feedback 
even in massive galaxies - otherwise, even in the absence of cool- 
ing, the gas coming from stellar recycling and mergers is enough 
to keep up relatively high star formation rates in massive galaxies. 



be higher than 2.5, and we enforce a constant value for /d 
at z > 5.7. The value of /d as a function of halo mass and 
time is shown in Fig. [S] top panel. 

• There is no reincorporation, i.e. we assume that cold 
gas that has been made unavailable for star formation by 
feedback once never returns back to the cold gas reservoir. 

• We treat dynamical friction as in NWIO, with a prefac- 
tor of Odf =3. 

• No merger-induced star bursts are included. 

• The gas residing in the filaments of satellites is stripped 
with an exponential timescale of 4 Gyr. 

Overall our model has similar basic scalings as other SAMs, 
although the parameter values are slightly different, and the 
recipes are simplified. 

We ran the standard model (as well as all following 
models) on the full volume of the Millennium Simulation 
(Springel et al. 2005). The resulting stellar mass functions 
(SMFs) at different redshifts and the sSFR-z relation (de- 
fined in what follows as the sSFR as a function of redshift 
at stellar mass 2 • 10^ - 10^'' Mo) are shown in Fig. [7] and 
[8] as thin blue lines. To account for the observational com- 
pleteness limit as indicated by Stark et al. (2009), we only 
take into account galaxies with log(SFR) > 0.25, 0.45 and 
0.5 yr~^ at z—4:, 5 and 6 respectively in Fig. [8] Clearly, the 
observed relation in the mass bin 2 • 10^ - 10^° M© is not 



© 2010 RAS, MNRAS 000, 000-000 



8 S. M. Weinmann, E. Neistein, A. Dekel 



reproduced, very similar to the models from NWIO shown 
in Fig. [2l On the other hand, note that the the evolution of 
the SMF is reproduced well (FigO, which is important for 
the following discussion of non-standard models. We show 
the relation between sSFR and stellar mass at z=4 and z=6 
in Fig. [9l top left panels. The evolution of the mean cos- 
mological SFR density is shown as thin blue line in Fig. 1101 
We have checked that including merger-induced bursts ac- 
cording to Croton et al. (2006) does not have a significant 
impact on any of the results shown here. 

4.2 Tuning 

We next gradually modify selected ingredients of the stan- 
dard model. Our first goal is to bring the median sSFR at 
2 = 2 — 6 and in the mass bin (0.2 - 1) • IO^V/q into bet- 
ter agreement with the observed plateau, the grey band in 
Fig. [T] We also consider the relation between sSFR and stel- 
lar mass at z = 4 and z = 6, shown in Fig. [21 in order to 
verify that the median sSFR tends to decrease with increas- 
ing stellar mass. We do not attempt a fit with the observed 
sSFR-mass relation outside the mass bin 
because the lower mass bin is strongly affected by incom- 
pleteness (e.g. Stringer et al. 2011), while results for the 
higher mass bin seem controversial (compare Stark et al. 
2009 to Lee et al. 2010a). Our secondary goal is to avoid se- 
vere deviations from the observed SMF at intermediate red- 
shifts, and from the observed evolution of mean SFR density. 
We do not attempt to fit the observed sSFR at z < 2. 

This is because we believe the mismatch in the sSFR 
between observations and models at 2; < 2, where the SFR 
in the universe declines rapidly, is a separate problem that 
deserves a dedicated study. It is discussed in more detail 
elsewhere (e.g. Fontanot et al. 2009; Karim et al. 2011). 
Additionally, we find fitting the z > 2 sSFR challenging 
enough even prior to adding the additional constraints im- 
posed by the low redshift data. On a separate note, we 
point out that observations of sSFR in galaxies with masses 
(0.2—1) • 10^° at 2 = 1 — 2 are rather uncertain, probably 
more so than at higher redshifts, as they tend to fall below 
the current stellar mass completeness limits at these red- 
shifts (e.g. Rodighiero et al. 2010, Dunne et al. 2009). This 
is due to the decline of the global star formation rate density 
and the increasing importance of dust, which decrease the 
intrinsic luminosity of galaxies at z < 2. 

For tuning the models, we focus on one model ingredient 
at the time, i.e. the star formation efficiency or the feedback 
efficiency. Using a table with discrete values of the efficiency 
in 8 bins in time, and 10 bins in halo mass, we start by fitting 
the plateau at the highest redshift and then subsequently 
continue to lower redshifts. We try to keep dependencies on 
halo mass and time as monotonic as possible, to limit the 
number of possible models, and also because monotonic de- 
pendencies are easier to motivate physically. Once we have 
found a solution for the plateau from 2; = 2 — 7, we compare 
the model results with the other observational constraints. 
Depending on the outcome, we discard the model or improve 
the fit to the stellar mass functions and the sSFR-stellar 
mass relation by making additional changes to the efficien- 
cies. Part of our models (SF2, SF3, FB2) are then further 
improved by simple changes to one or two other model ingre- 
dients. We point out that we have not carried out a system- 



Table 1. A summary of the models discussed in this work. Eb is 
the star formation efficiency, the feedback efficiency, /rc the 
reincorporation efficiency, /burst the efficiency of merger-induced 
star bursts, and a^f is the dynamical friction prefactor. mi and 
m2 is the sum of the stellar and cold gas mass in the merger 
progenitors. If not listed, elements are kept at the standard model 
values. 



Models Modifications line type 



standard - thin blue 

SFl es(t) thin dashed red 

SF2 fs(i), /burst (ti ™i/™2) dotted-dashed green 

SF3 es(t), /burst (t, "11 /"i-2), a!df(t) dotted brown 

SF4 es(t, Miiaio) dashed cyan 

FBI /d(t, A^halo) dashed pink 

FB2 /d(t,A4alo), /rc(t,A4alo) thick black 



atic study covering all the parameter space, but a process 
of trial and error that is geared towards finding successful 
solutions which fit the constraints in question. Fitting the 
sSFR plateau and the evolution of the stellar mass function 
together is not trivial, and a substantial number of itera- 
tions were needed for each model to arrive at the solutions 
presented below. 

The different models are summarized in table [T] 

4.3 Reproducing the sSFR plateau 

In this section, we explore two alternative modifications 
to the standard model that aim at reproducing the sSFR 
plateau: (i) models "SF" , in which Cs is reduced at high red- 
shifts, and (ii) models "FB", where the feedback efficiency 
is enhanced at high redshifts. Results for all the models are 
shown in Fig. [7l[8l[9l andfTOl 

4.3.1 Model SFl - tune ts 

In model SFl, we tune the SFE parameter, es, while all other 
model parameters are kept fixed as in the standard model. 
With 65 varying in time as shown in Fig. [S] (top panel) , the 
model reproduces the sSFR plateau at 2 = 2 — 6 as shown 
in Fig. [8] While es does not need to depend on halo mass, 
it is not monotonic with time. In order to have an sSFR 
plateau starting from redshift Zp (chosen here to be 2p ~ 7), 
this model requires an earlier epoch where ts is well above 
its fiducial value 0.01. This early star formation is needed 
in order to produce enough galaxies of large-enough stel- 
lar mass by Zp, after which the much-lower SFR adds only 
slowly to the stellar mass. At the onset of the plateau, has 
to drop to values well below 0.01, and it should continue to 
gradually decline till 2 ~ 4, in order to permit the observed 
low sSFR values at 4 < 2 < 2p. After 2 ~ 4, ts is gradually 
increasing in order to match the high observed sSFR in this 
regime. It catches up with the fiducial value es ~ 0.01 only 
at low redshifts. We note that model SFl predicts that the 
sSFR plateau does not extend all the way to the epoch of the 
emergence of the first galaxies; it is preceded by a starburst 
epoch. 

Unfortunately, Fig. [7] shows that model SFl does not 
produce enough massive galaxies to match the observed 



© 2010 RAS, MNRAS 000, 000-000 



9 




Figure 7. Stellar mass function (SMF) of galaxies at different redshift bins, for our different models as indicated and as summarized in 
tablc[l] Observational data are marked by grey symbols. At 2 > 0, the model stellar mass is convolved with a Gaussian error distribution 
of standard deviation 0.25 dex, which includes the differences in the IMFs assumed by the different observers. Data at z = are by Li et 
al. (2009, circles), Baldry et al. (2008, crosses), and Panter et al. (2007, pluses). Data at z > are from Bundy et al. (2006, z = 0.75 — f , 
circles), Borch et al. (2006, z = 0.8 — f , crosses), Perez- Gonzalez et al. (2008, z = 0.8 — I, z = 1.6 — 2, 2 = 2.5 — 3, plus signs), Fontana 
et al. (2006,2 = 0.8 - 1, z = 1.6 - 2, 2 = 2 - 3, stars), Drory et al. (2004, 2 = 0.8 - 1, upward-pointing triangles), Drory et al. (2005, 
z = 0.75 — 1.25, 2 = 1.75 — 2.25, 2 = 2.25 — 3 , diamonds and squares), Marchesini et al. (2009, 2 = 1.3 — 2, 2 = 2 — 3, right-pointing 
triangles). Model SMFs are plotted at 2 = (top left panel), 2 = 0.8, 1, 1.2 (top right panel), 2 = 1.2, 1.5, 2 (bottom left panel) and 2=2, 
2.5, 3 (bottom right panel). 




Figure 8. Evolution of sSFR for galaxies with stellar mass in the fixed bin 2 ■ 10 — 10 Mq for the models indicated and listed in table 
[1] At 2 > 4, only galaxies with SFR above the completeness limits given by Stark et al. (2009) are included. Errorbars denote ±lcr at 
2 > 2. Each of the redshift bins contains more than 50 galaxies. All the models shown reproduce the sSFR plateau at 2 = 2 — 6. 
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Figure 9. sSFR as a function of stellar mass, at z=4 (top panels) and z = 6 (bottom panels). The results of model FBI are very similar 
to model FB2. The greyscale refers to log number density of galaxies in units of Mpc~^dex~^. 



mass function at intermediate redshifts, despite the high ini- 
tial Es. Model SFl thus reveals a nontrivial tension between 
two sets of data, namely the low values of sSFR at high red- 
shifts versus the high-mass end of the SMF at intermediate 
redshifts. This indicates that the high values of sSFR ob- 
tained in SAMs at high redshifts are needed there for the 
purpose of building up massive enough galaxies by 2 ~ 2, 
and are therefore not easy to avoid. 



4.3.2 Model FBI - tune feedback 

In model FBI, we only tune the stellar feedback parameter 
fd , while all other model parameters are kept fixed as in the 
standard model. The values of fd as a function of mass and 
redshift are shown in Fig. |6] (top-right panel) , and the resul- 
tant sSFR evolution and SMF are shown in Fig.[7]and[8l We 



find that fitting the sSFR plateau requires very high stellar 
feedback efficiencies at early times (z >, 3), and especially 
so for massive haloes. This is needed in order to balance 
the high accretion rates of these haloes. High feedback ef- 
ficiency is not needed in the massive haloes at later times 
because the cooling in them is set to zero at 2 > 3. At lower 
halo masses Mh ^ IO'^^Mq, moderately strong feedback is 
needed at all redshifts. Note that stellar feedback is not en- 
tirely monotonic here with either mass or time; we improve 
this in the following model FB2 by adding reincorporation 
as an additional model ingredient. 

As in model SFl, we do not reproduce the high mass end 
of the SMF with this model, again indicating that suppress- 
ing the efficiency of star formation at high redshift leads to 
an underproduction of massive galaxies, if no other changes 
to the model are made. As can be seen in Fig. [T] this short- 
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Figure 10. Cosmological evolution of SFR density for mod- 
els as indicated versus observations. The upper grey belt repre- 
sents the compilation by Hopkins & Beacom (2006, Scr confidence 
level) from measurements of SFR. More recent determination by 
Bouwens et al. (2009) and by Kistler et al. (2009) from gamma-ray 
bursts are consistent with the upper grey belt. The lower grey belt 
is from Wilkins et al. (2008; 1(7 confidence level), derived from the 
growth curve of the stellar mass density. The model predictions 
lie in between. 

age is more severe than in the SF models, and more so at 
higher redshifts. We did not manage to improve this aspect 
of the FBI model so far, since tuning the FB models is more 
difficult than tuning the SF models. The reason for this is 
that the star formation rate of galaxies is directly propor- 
tional to es, while its dependence on the feedback efficiency 
is non-linear. 

All our models that reproduce the sSFR plateau at 
z < 7 require an earlier phase of high SFR, with a sSFR 
higher than the plateau level by a factor of a few. This pro- 
vides a high enough stellar mass at z ~ 6, which allows the 
desired low sSFR at the plateau level while the SFR is driven 
to high values by the high accretion rate. For example, in 
order to reach by z ~ 7 a mass of 2 • 10^ Mq, the minimum 
mass of galaxies on the plateau, a galaxy needs an average 
SFR ~ 10 M0yr~'^ between the onset of star formation at 
z ~ 9 and z ~ 7, say. This implies a sSFR ~ 5Gyr~^ at 
~ 7, and even higher values at earlier times, when the 
stellar mass is smaller. Such an early phase of high SFR has 
to be introduced by hand when es is set to low values at 
z < 7, as in model SFl. In model FBI, on the other hand, 
the high SFR at 2 > 7 occurs naturally, because the gas 
available for SFR in a growing galaxy is proportional to the 
instantaneous accretion rate, which is only slowly increasing 
with time. This is in contrast to the SF models, where the 
SFR is proportional to the accumulated gas in the galaxy 
and therefore tends to increase faster with time. 



4-3.3 Model FB2 - tune feedback and reincorporation 

In model FB2, we improve on model FBI by adding rein- 
corporation from the feedback phase back to the cold phase 
as an additional model ingredient. This makes it possible to 
keep the variation of the feedback efficiency with time mono- 
tonic (although the trend is of opposite sense at low and high 
halo masses). With model FB2, we can reproduce the sSFR 



plateau very well down to 2: --^ 2 (Fig. [S] black thick solid 
line), but again fail to simultaneously reproduce the high 
mass end of the SMF at z >, 2 (Fig.[71 black thick solid line). 
Feedback and reincorporation efficiencies in model FB2 are 
shown in Fig. (6] left and right bottom panel respectively. 
Reincorporation is needed in order to keep up high sSFR 
at « < 3 for haloes with masses Afh ^ 10^^ M©. We set it 
to very low values outside this range. The efficiency of rein- 
corporation in the range of halo mass and time when it is 
needed then quite similar at all redshifts when expressed in 
units of 1/tdyn, with a slight peak at around z~ 2. 



4.4 Reproducing the sSFR and the evolution of 
the SMF 

While we have found in the previous section that it is possi- 
ble to reproduce the sSFR plateau once the star-formation 
or feedback efficiency is allowed to vary in time and with 
mass, all the models discussed so far underproduce the high 
mass end of the SMF at z >, 2. In this section, we attempt 
to improve the fit to the SMF at intermediate redshifts by 
an additional modification to the model. For simplicity, we 
focus on modifications to model SFl. Figure [7] shows that 
these models, SF2-SF4, reproduce the high mass end of the 
intermediate-redshift SMF better than the previous models. 
Figure [§] demonstrates that this success is associated with a 
population of galaxies more massive than l(f Mq with high 
SFR at z ~ 4, which were missing in the previous mod- 
els. We explore here three different ways for producing this 
population. In model SF2, we boost the efficiency of merger- 
induced star bursts. In model SF3, we speed up the merger 
rate. In model SF4, we introduce a rather involved variation 
of both with halo mass and time, making star formation 
more efficient in relatively high mass haloes while keeping it 
inefficient in lower mass haloes. 



4-4- 1 Model SF2 - tune and mergers 

In model SF2, we reproduce the required massive galaxies 
at 2: ~ 2 by boosting up merger-induced starbursts. We first 
tried the starburst prescription of Croton et al. (2006) where 
/burst = 0.56(7711/7712)" '^, but this had only a little effect on 
the SMF. However, with /burst ~ (7771/7772)" '' at z > 1, 
followed by the Croton et al. (2006) prescription at 2 < 
1, model SF2 produces a higher abundance of high mass 
galaxies in better agreement with the SMF at intermediate 
redshift (Fig. [7]), while still reproducing the sSFR plateau 
(Fig.®. 

An additional parameter governing merger- induced 
starbursts is the burst duration, which is set to 10 Myr. 
We find that smearing the bursts over a much longer du- 
ration (e.g. 500 Myr) makes very little difference to the re- 
sults. As we only consider the median and the 68% range, 
a small population of starbursting galaxies can be present 
without ruining the sSFR plateau. Strongly star bursting 
galaxies might also not make it into the Stark et al. (2009) 
or Gonzalez et al. (2010) samples, as the highly star form- 
ing regions in galaxies tend to be strongly dust attenuated 
(e.g Chariot & Fall 2000), and are therefore a good way of 
building up stellar mass in a hidden mode. 
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Model SF3 - tune eg and mergers 

Model SF3 is a variation of model SF2, where the high star- 
burst efficiency at 2 > f is replaced by a shorter charac- 
teristic time for dynamical friction. The fiducial dynamical 
friction prefactor of Qdf = 3 is replaced by Qdf = 0.1 at z > 1 
and ftdf = 5 at 2 < 1, while the starbursts are moderate, 
/burst ~ 0.2(mi/m2)'''^. The rapid stellar assembly in this 
model boosts up the buildup of stellar mass, and we need 
to lower the efficiency of star formation accordingly (Fig. O 
top panel). This model provides a sensible fit to the SMF 
(Fig. [7|) and it reproduces the sSFR plateau (Fig. |8]). 

4.4.3 Model SF4 - tune 

Model SF4 demonstrates a third way to build up high mass 
galaxies by moderate redshifts without enhancing the con- 
tribution of mergers. In this model, the star-formation ef- 
ficiency, es, is varied as a function of both halo mass and 
time as shown in Fig. [5] This allows for a simultaneous fit 
to the SMF (Fig. [3 and the sSFR plateau (Fig. [8]). We note 
that in this model, the sSFR does not in general decrease 
with increasing stellar mass, as shown in Fig. |9] for the z=6 
results, bottom right panel. This is contrary to all obser- 
vations of this relation we are aware of, which tend to find 
that the median sSFR of galaxies decreases with increasing 
stellar mass. 

4.5 Summary of model results 

In this section, we summarize the successes and failures of 
the models, guided by Fig. [7l|8]and [10] From Fig. [7] which 
shows the stellar mass functions up to 2 ~ 3, it becomes 
clear that only models SF3 and SF4 manage to reproduce 
the SMF at z> 1.2. 

In Fig. Is] we show that all our models except the stan- 
dard models provide a reasonable fit to the sSFR plateau 
from 2; ~ 2 — 6. At z < 2, the sSFR is lower than the ob- 
servational estimates. This may be partially due to incom- 
pleteness in the observations, which misses galaxies with low 
SFR especially at low redshift, where the global SF in the 
Universe has declined. The mismatch at z ~ 0, however, is 
probably real and already apparent in the standard model. It 
is likely connected to a similar underproduction of the sSFR 
present in other recent SAMs (e.g. Fontanot et al. 2009, Guo 
et al. 2011). We note that our efforts to improve the match 
to the sSFR plateau lead to an increase in the sSFR at 2 < 2 
and thus to a better agreement with observations. 

In Fig.[lO]we show the global star formation rate density 
as a function redshift in our models. The upper grey belt 
represents the directly observed SFR density from Hopkins 
& Beacom (2006), while the lower belt is derived from the 
growth curve of the stellar mass (Wilkins et al. 2008). More 
recent estimates like from Bouwens et al. (2009) and from 
Kistler et al. (2009) are in agreement with the upper belt. 
The observational results that are directly measured from 
star formation, and those that are obtained indirectly from 
the evolution of the stellar mass density are thus clearly 
inconsistent, which is the reason that SAMs (like e.g. Guo 
et al. 2011) usually have a SF density somewhere in between 
these two constraints. Results of our various models start to 
deviate considerably at 2; > 2, but the discrepancy between 



the observational results makes it impossible to use them to 
constrain models. The highest global SFR density at 2 > 3 
is reached by the standard model, followed by the models 
which have the best fit to the SMF at high redshift, namely 
model SF3 and SF4. 



4.6 The sSFR of individual galaxies 

In Fig. 111! we follow the growth of one individual model 
galaxy following the main branch of its merger history. 
Shown as a function of redshift is the mass in the dark- 
matter, cold gas and stars, the SFR and sSFR, for the stan- 
dard model, SF2 and FB2. The resulting stellar mass at 
2 = is similar in the three models. However, the star- 
formation history and the growth of stellar mass are differ- 
ent — in both model SF2 and FB2, the buildup of stellar 
mass is delayed compared to the standard model. 

We note that at 2 > 2 the SFR in individual galaxies 
is rising with time for all models. This is a common feature 
to all our model galaxies, and is consistent with the pre- 
dictions of Finlator et al. (2006, 2011), as well as with the 
observational finding of Papovich et al. (2010), Maraston et 
al. (2010), Lee et al. (2010a,b). It is emerging naturally from 
the fact that in an individual galaxy, the accretion rate is 
expected to be growing in time in its rapid growth phase at 
2 > 2 (Neistein & Dekel 2008). 

We see in Fig. [11] that as the galaxy evolves from 2 = 6 
to 2 = 2, the sSFR declines in a rather slow pace under 
models SF2 and FB2 compared to the steeper decline in the 
standard model. The sSFR plateau discussed in the previous 
sections, which refers to the median value for galaxies of a 
fixed mass at different redshifts, does not necessarily imply 
a constant sSFR in individual galaxies as they grow. How- 
ever, the observations of Stark et al. (2009) and Papovich 
et al. (2010, Figs. 2 and 3) do argue for a constant sSFR in 
individual galaxies as well. It is thus encouraging that our 
models that reproduce the sSFR plateau also come close to 
reproducing the constancy of the sSFR in individual galax- 
ies. 



5 DISCUSSION 

In subsections l5.ll and l5.2l of this discussion section, we com- 
ment on the physical plausibility of the ad-hoc modifications 
suggested by our models. This discussion is not complete or 
conclusive — it is only meant to raise some of the relevant 
issues and to trigger further discussion. The limited goal of 
the current paper remains to point out the possible nature of 
the modifications to the standard models that are required 
for reproducing the sSFR plateau together with the SMF at 
moderate redshift. 



5.1 Changes to the quiescent evolution of galaxies 

(i) Star formation efficiency: In the SF models, Cs has to 
vary in time, with a very high efficiency at 2 > 7, and a low 
efficiency at 2; = 3 — 6. The latter is consistent with the pre- 
diction of Gnedin & Kravtsov (2010) that the low metallicity 
and high UV radiation at 2 ~ 3 should significantly lower the 
normalization of the Kennicutt-Schmidt relation. This effect 
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Figure 11. Evolution in a single galaxy within the same halo according to different models. Top row: standard model. Middle row: 
model SF2. Bottom row: model FB2. Left: mass in dark matter (same in all models), stellar mass, and cold gas. Right: SFR and sSFR. 
Horizontal dashed lines indicate the stellar mass range relevant for the median sSFR shown in Fig. [S] The peaks in the SF history of 
model SF2 are due to starbursts. 



is particularly strong for low gas surface densities. An effec- 
tive suppression of star formation in high redshift galaxies 
and at low masses turns out to be a natural outcome of the 
low metallicity in these galaxies (Krumholz & Dekel 2011). 
The low dust content enables the UV radiation from stars to 
heat the gas (while dissociating H2 molecules) and to pre- 
vent further star formation, until enough metallicity is built 
up and the SFR regions are shielded. Observational studies 
at z ~ 3 seem to agree with this prediction (e.g. Wolfe & 
Chen 2006; Rafelski et al. 2010). Furthermore, both SAMs 
and hydrodynamical simulations have suggested that a low 
Eg at high redshift improves the match with observations in 
a number of ways (Baugh et al. 2005; Lacey et al. 2010; 
Agertz et al. 2011). It is more difficult to explain the high 
Es at z > 7, necessary for producing enough galaxies with a 
high stellar mass by 2 ~ 6. What might help is a feedback 



mechanism that causes an abrupt change in the mode of 
star formation after a very active initial phase. Model SF4 
suggests in addition an increase of with halo mass. Such a 
variation might result from faster metal enrichment or less 
efficient feedback in more massive haloes. 

(ii) Stellar feedback: The FB models suggest strong early 
stellar feedback that is not limited to small galaxies. This 
is in potential disagreement with conventional models of SN 
feedback, in which the effect of feedback at a fixed halo 
mass is expected to be weaker at high redshift, where the 
deeper potential well makes it harder to eject gas from the 
galaxy (Dekel & Silk 1986; Dekel & Birnboim 2006). One 
should therefore consider additional feedback mechanisms 
that may vary in the required way with time and mass. For 
example, a metallicity-dependent stellar feedback (Nishi & 
Tashiro 2000; Krumholz & Dekel 2011) will be especially 
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effective at high redshifts. The required strong feedback at 
high halo masses (log(Mh) > 11.5 at tnubbic < 2) may be 
achieved by AGN feedback from massive black holes rather 
than stellar feedback. 

(iii) Reincorporation: Model FB2 requires a high rate of 
reincorporation only at z <, 3. In other SAMs, the reincor- 
poration rates scale with l/t^yn, and therefore gradually 
increase towards higher redshift, making the effective stel- 
lar feedback less efficient at high z. A weaker time depen- 
dence, more consistent with FB2, is found in simulations 
that incorporate strong momentum-driven winds (Oppen- 
heimer et al. 2010). Another option is "reincorporation" that 
utilizes a gas reservoir in the galaxy rather than gas that has 
been ejected earlier. For example, this occurs naturally for 
met allicity- dependent feedback. As long as the metal con- 
tent in a galaxy is low, gas is not well shielded from UV 
radiation by dust, leading to efficient heating and H2 disso- 
ciation. Once the metal content reaches a certain threshold 
value, the gas is shielded and star formation becomes effi- 
cient. The abrupt appearance at z ~ 3 of reincorporation 
in FB2 is an oversimplification, but the general increase of 
reincorporation with halo mass is plausible both for ejective 
and non-ejective feedback. 



5.2 Changes to the recipes for mergers 

(i) Enhanced starburst efficiency at high redshift: In 
model SF2, the starburst efficiency, defined according to eq. 
1131 is higher at higher z. This is not to be confused with the 
SF efficiency Cs, which is relevant for quiescent star forma- 
tion. An increased star burst efficiency at high z is plausible 
due to the shorter dynamical times. The gas-rich mergers 
at high redshift may in principle be very different from to- 
day's more familiar gas-poor mergers. High-resolution gas- 
rich merger simulations indeed indicate an effective star for- 
mation efficiency that is higher than in non-merger situa- 
tions (Teyssier, Chapon & Bournaud 2010). 

(ii) Enhanced merger rate at high redshift: In model SF3, 
the dynamical friction prefactor Odt is reduced to 0.1 at 
z > 1, which speeds up the merger process by a factor 
~ 20 — 30 compared to the standard model and other SAMs 
(Croton et al. 2006; De Lucia et al. 2007). At high redshift, 
the merging galaxies flow from the virial radius to the cen- 
tral galajcy along narrow radial streams associated with the 
cosmic-web filaments in about one half of a halo crossing 
time, -Rvir/Vvir, which is ~ O.ltuubbiG (Dekel et al. 2009). 
This implies a merging time substantially shorter then the 
dynamical friction time estimated for gradual spiraling in. 
Hopkins et al. (2010) point out that the standard estimate of 
merger time based on dynamical friction is indeed an over- 
estimate when the approach of the satellite is along a radial 
orbit, or when the dynamical friction estimate starts at a 
distance larger than 0.1 - 0.2 -Rvir, where the satellite is no 
longer properly resolved in the simulation. Both of these 
conditions tend to be valid at high redshift (Wetzel 2010, 
Hopkins et al. 2010), and together they may lead to an or- 
der of magnitude overestimate of the merger time in the 
standard model. We used the Millennium simulation to ver- 
ify that with Odf = 0.1, the merging time is comparable to 
the halo crossing time. 



5.3 Dust at high redshift 

The validity of the sSFR plateau at high redshifts crucially 
depends on the correction adopted for dust extinction by 
the different authors. At z ~ 2, the estimates for dust ex- 
tinction in LBGs from the UV slopes are confirmed by ra- 
dio estimates (Pannella et al. 2009) and from comparison 
with a local sample of Lyman break analogs (Overzier et al. 
(2011). At higher redshifts, the estimates are more uncer- 
tain. Stark et al. (2009), Gonzalez et al. (2010) and Labbe 
et al. (2010a,b) all assume practically no dust extinction at 
2: > 4. This assumption is supported by simple theoretical 
considerations and several observational studies. It is ex- 
pected that high-redshift galaxies should have a lower dust 
content than today's galaxies simply because they had less 
time to produce metals and dust. Indeed, theoretical mod- 
els (e.g. Guo & White 2009) often assume a lower dust ex- 
tinction at higher redshift. Bouwens et al. (2009) estimated 
the dust extinction at high redshift based on the observed 
UV continuum slope, and found low values at 2: >, 4, espe- 
cially for the dominant population of galaxies with relatively 
low UV luminosities. Labbe et al. (2010b), Finkelstein et 
al. (2010) and Bouwens et al. (2010) also report on low or 
zero dust extinction at z ~ 6 — 8. Brammer & van Dokkum 
(2007) tested whether the usual LBG selection criteria miss 
a significant number of dusty galaxies at z ~ 4 by selecting 
according to Balmer breaks instead, and found that this is 
not the case, again indicating that strong dust obscuration 
at high redshift is rare. This may however not be the case for 
massive galaxies, which potentially already have high dust 
content at 2 ~ 3 — 4 (Mancini et al. 2009, Marchesini et 
al. 2011) and thus may be missed by the usual LBG selec- 
tion criteria. Also, Schaerer & de Barros (2010) and Yabe 
et al. (2009) argue for high correction factors of ~ 9 when 
translating the UV fiux to SFR at 2 = 5 — 7. We note that 
in order to reconcile the plateau in the sSFR with current 
SAMs by the effect of dust-extinction alone, the dust extinc- 
tion would have to increase with redshift, a trend which is 
opposite to basic theoretical expectations. While the "dust" 
has not settled yet on this debate, it seems that the evidence 
for a sSFR plateau is intriguing enough to justify a serious 
theoretical consideration, but it is left for future observa- 
tions to tell whether this is indeed a valid strong constraint 
or a fiuke. 



5.4 Comparison to previous work 

Previous efforts to understand the properties of star forming 
galaxies at 2 > 4 include Finlator et al. (2006, 2011), Night 
et al. (2006), Mao et al. (2007), Nagamine et al. (2008), 
Stark et al. (2009), Lee et al. (2009), Khochfar & Silk (2011), 
Lacey et al. (2010), Lo Faro et al. (2009), and Stringer et 
al. (2010). Here we compare our findings to some of those 
studies. 

Lee et al. (2009) concluded based on clustering analysis 
that the duty cycle of star formation in high-redshift galaxies 
must be short, consisting of episodes of < 0.35 Gyr in each 
galaxy. Stark et al. (2009) came to similar conclusions. This 
would be in apparent conflict with the need to build up the 
massive end of the SMF by 2 ~ 2 - 3. Our model SFl that 
reproduces the sSFR plateau assumes a continuous mode of 
star formation, and even then, it underpredicts the massive 
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end of the SMF. If the low sSFR were associated instead with 
short-Uved bursts of star formation, it would have been even 
harder to build up the stellar mass quickly enough. A way to 
reconcile a short duty cycle with the required fast build-up 
of stellar mass might be that the short-lived episodes of SF 
are only the final phase of a much more active and dust- 
obscured starbursts, enough for building up the high mass 
end of the SMF but undetectable in their violent phase. 
Finlator ct al. (2011) suggest that the clustering data by 
Lee et al. (2009) can also be explained without a bursty star 
formation history, if there are strong outflows that lead to 
an increased scatter in the relation between baryonic mass 
and host halo mass. 

Khochfar & Silk (2011) present a model that reproduces 
the sSFR plateau at « > 4 with (i) a very high frequency of 
star-bursts at high redshift and (ii) a burst efficiency that 
scales with the inverse of the halo circular velocity to the 
third power. In this way, star formation in galax;ies in the 
mass range relevant for the plateau is inefficient at high red- 
shift, but stellar mass grows efficiently during mergers in low 
mass haloes. Their model thus resembles a combination of 
our models SF2 and SF3. 

The simultaneous need for a low sSFR at high redshift 
and a fast build-up of stellar mass seems to be in odds with 
the notion that the observed growth rate in stellar mass 
density is low compared to that predicted by integrating the 
observed SFR density (Wilkins ct al. 2008). These seemingly 
inconsistent problems may reflect a difference between the 
overall evolution of stellar mass density and the evolution of 
the population dealt with here, involving relatively massive 
galaxies in a fixed mass bin. 

Lo Faro et al. (2009), comparing a SAM for LBGs at 
z ~ 4 — 6 to observations, found that their model pro- 
duces an excess of star-forming galaxies with stellar masses 
~ 10^ - 10^" Mq at 2 ~ 4, within the range relevant for 
our study of the sSFR plateau. They conclude that some 
feedback mechanism must suppress star formation at early 
times in the corresponding haloes of ~ 10^^ — lO^^M©, but it 
should become ineffective at later times in haloes of similar 
mass. As Lo Faro et al. (2009) point out, such a behaviour of 
the feedback efficiency is not part of the standard feedback 
models. Their conclusions are thus similar to ours, albeit 
they originate from a different observational constraint. 



6 CONCLUSIONS 

We explored possible ad-hoc modifications to standard semi- 
analytic models that reproduce a constant sSFR at 2; > 2 
while growing enough massive galaxies by ~ 2, and came 
to the following conclusions. 

• We confirm that a sSFR plateau ai z — 2 — 6 is in robust 
disagreement with current models, as claimed in Bouchc ct 
al. (2010). Wc have demonstrated that, in a fixed stellar 
mass bin of 2 ■ 10^ — 10^" Mq , the common feature of a variety 
of standard SAMs is a gradual decline of the sSFR with time, 
associated with the decline of the total specific accretion 
rate. 

• We find that it is possible to reproduce the sSFR 
plateau together with the stellar mass function at ~ 2 
via non-trivial modifications to the standard SAMs. Three 



different modifications seem necessary, related to three dif- 
ferent observational features: (1) the low sSFR at 2 > 4, (2) 
the high sSFR at « ~ 2 — 3, and (3) the abundance of high 
mass galaxies at 2 2 — 3. The low sSFR at 2 > 4 is repro- 
duced either by strong stellar feedback at high redshift in 
all masses, or by inefficient star formation at high redshift 
following a phase of very efficient star formation at very high 
redshift. The high sSFR at 2 = 2 — 3 could emerge either 
by a drop in the feedback efficiency at this epoch, or by a 
corresponding enhancement of star formation efficiency, or 
by efficient reincorporation of gas that was previously pre- 
vented from forming stars. Finally, the high mass end of the 
SMF at 2 > 1 can be generated despite the low SFR at high 
redshift by an additional modification of the star formation 
in a sub-population of massive galaxies at high redshift. This 
is achieved in our models by either speeding up the mergers, 
or enhancing the merger-induced starbursts, or by introduc- 
ing a non-trivial dependence of star formation efficiency on 
halo mass. However, none of our modified models matches 
the stellar mass function at 2 ~ 3 as well as our standard 
model. This refiects the fact that the low sSFR at high red- 
shift and the presence of massive galaxies by 2 ~ 3 are not 
easily reconciled. 

• Our models predict that the SFR in individual galax- 
ies is monotonically increasing with time at 2 > 2. This is 
in agreement with the theoretical predictions of Finlator et 
al. (2006, 2011) and Bouche et al. (2010), as weU as the 
observational results by Papovich et al. (2010). We note in 
particular that an SFR that grows exponentially with time 
implies a constant sSFR. It should be mentioned that the 
decreasing SFR sometimes assumed in SED fitting (e.g. in 
Stark et al. 2009) is incorrect. 



Wc have demonstrated that the simple SAM of NWIO is 
useful for exploring how current SAMs should be modified in 
order to match new observational constraints, and for point- 
ing out apparently conflicting observational constraints. A 
similar method will be useful in addressing other puzzling 
observations such as the tilt in the relation of sSFR and stel- 
lar mass (e.g. Somerville et al. 2008), the fraction of passive 
galaxies as a function of stellar mass at 2 = (e.g. Wein- 
mann ct al. 2010), or the fraction of AGNs as a function of 
stellar mass (Fontanot et al. 2010). 

We learn that the observed sSFR at high redshift has 
the potential for posing powerful constraints on the phys- 
ical processes of star formation and feedback. If the sSFR 
plateau, as observed by Stark et al. (2009), Labbe ct al. 
(2010a, 2010b), and Gonzalez et al. (2010), is confirmed, it 
will provide invaluable information on the baryonic physical 
processes at high 2, indicating that it could be different from 
those at low redshift. One should be eagerly waiting for new 
developments in the observations of SFR and stellar mass 
at high redshift. 

While we started a discussion of the physical plausibil- 
ity of the required modifications to the standard recipes of 
galaxy-formation models, the modifications suggested above 
are primarily of an ad-hoc nature. They deserve a thorough 
theoretical study of physical mechanisms that could be re- 
sponsible for the required variation with time and mass. 
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